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Jets and environment of microquasars 



S. N. Zhang 1 ; J. F. Hao 1 

Abstract. Two relativistic X-ray jets have been detected with the Chandra 
X-ray observatory from the black hole X-ray transient XTE J1550-564. We 
report a full analysis of the evolution of the two jets with a gamma-ray burst 
external shock model. A plausible scenario suggests a cavity outside the central 
source and the jets first travelled with constant velocity and then are slowed 
down by the interactions between the jets and the interstellar medium (ISM). 
The best fitted radius of the cavity is ~0.36 pc on the eastern side and ~0.46 pc 
on the western side, and the densities also show asymmetry, of ~0.015 cm~ 3 on 
the east to ~0.21 cm -3 on the west. A large scale low density region is also found 
in another microquasar system, H 1743-322. These results are consistent with 
previous suggestions that the environment of microquasars should be rather 
vacuous, compared to the normal Galactic environment. A generic scenario 
for microquasar jets is proposed, classifying the observed jets into three main 
categories, with different jet morphologies (and sizes) corresponding to different 
scales of vacuous environments surrounding them. 



1. INTRODUCTION 

Microquasars are well known miniatures of quasars, with a central black hole 
(BH), an accretion disk and two relativistic jets very similar to those found in 
the centers of active galaxies, only on much smaller scales (Mirabel & Radn'guez 
1999). Since discovered in 1992, radio jets have been observed in several BH 
binary systems and some of them showed apparent superluminal features. In 
the two well known microqusars, GRS 1915+105 (Mirabel &; Radriguez 1999) 
and GRO J1655-40 (Tingay et al.1995; Hjellming & Rupen 1995), relativistic jets 
with actual velocities greater than 0.9c were observed. In some other systems, 
small-size "compact jets", e.g. Cyg X-l (Stirling et al. 2001), and large scale 
diffuse emission, e.g. SS433 (Dubner et al. 1998), were also detected. 

XTE J 1550-564 was discovered with RXTE during its strong X-ray outburst 
on September 7, 1998 (Smith 1998). It is believed to be an X-ray binary system 
at a distance of ~5.3 kpc, containing a black hole of 10.5±1.0 solar masses and 
a low mass companion star (Orosz et al. 2002). Soon after the discovery of the 
source, a jet ejection with an apparent velocity greater than 2c was reported 
(Hannikainen et al. 2001). In the period between 1998 and 2002, several other 
outbursts occurred but no similar radio and X-ray flares were detected again in 
these outbursts (Tomsick et al. 2003). 
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With the help of the Chandra satellite, Corbel et al (2002) found two large 
scale X-ray jets lying to the east and the west of the central source, which 
were also in good alignment with the central source. The eastern jet has been 
detected first in 2000 at a projected distance of ~2l" from the central black 
hole. Two years later, the jet could only be seen marginally in the X-ray image, 
while a western counterpart became visible at ~22" on the other side. The 
corresponding radio maps are consistent with the X-ray observations (Corbel et 
al. 2002). 

There are altogether eight two-dimentional imaging observations of XTE 
J 1550-564 in Chandra archive during June 2000 and October 2003 (henceforth 
observations 1~8). Here we report a full analysis of these X-ray data, together 
with the kinematic and spectral evolution fittings for all these observations. 

2. OBSERVATIONS of XTE J1550-564 

The basic information of observations 1~8 is listed in Table 1, including the 
observation ID, date, and the angular separation between the eastern and west- 
ern jets and the central source. The positions are obtained by the Chandra 
Interactive Analysis of Observations (CIAO) routine wavdetect (Freeman et al. 
2002). In observations 5 and 6, no X-ray source is detected by wavdetect at the 
position of the eastern jet. However, from the smoothed images (Fig.l), a weak 
source could be recognized in observation 6. We thus select the center of the 
strongest emission region as one data point. We calculate the source centroid for 
the central source and the X-ray jet respectively and for all the five observations, 
the calculated position changed by less than 0.5". Therefore, an upper limit of 
0.5" is set for the error of the jet distance. 

From Table 1 and Fig.l, we could see clearly that an X-ray emission source 
is detected to the east of the central source in the first four observations and 
another source is detected to the west in the last five observations. Calculations 
also show that these two sources, when presented in a single combined image, are 
in good alignment with the central compact object with an inclination angle of 
85.9°±0.3°. By calculating the average proper motion, an approximate estimate 
of deceleration could be seen for both jets. 

3. ENERGY SPECTRUM and FLUX 

Since the emission from the eastern jet has been studied fully (Corbel et al. 2002; 
Tomsick et al. 2003), we mainly focus our spectral analysis on the western jet. 
The X-ray spectrum in 0.3-8 keV energy band is extracted for each observation of 
the western jet. We use a circular source region with a radius of 4", an annular 
background region with an inner radius of 5" and an outer radius of 15", for 
each observation. Instrument response matrices (rmf) and weighted auxiliary 
response files (warf) are created using CIAO programs mkacisrmf and mkwarf, 
and then added to the spectra. We re-bin the spectra with 10 counts per bin 
and fit them in Xspec. 

The results of spectra fitting with an absorbed power-law model are also 
shown in table 1. We use the Cash statistic since it is a better method when 
counts are low. The absorption column density is fixed to the Galactic value 
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Figure 1. Left: The Chandra 0.3-8 keV raw image showing XTE J1550-564 
and the western jet. The images are in linear scale and no count saturation 
has been set. Right: The smoothed Chandra X-ray images of the eight ob- 
servations of XTE J1550-564 and the two jets together. The green elliptical 
regions are source emission regions by warded. Observation 4 shows the good 
alignment of the two jets and the central source. 



in the direction of XTE J1550-564 obtained by the radio observations (Njj = 
9 x 10 21 cm -2 ) (Dickey & Lockman 1990). Our results are quite consistent with 
previous work by Karret et al.(2003). The calculated absorbed energy flux in 
0.3-8 keV band is comparable to the value of the eastern jet. The observed flux 
decayed rather quickly, from ~ 1.9 x 10 _13 erg cm -2 s _1 in March 2002 to only 
one sixth of this value in October 2003 (see section 4.2). 



Table 1. XTE J1550-564 Chandra Observations 



Num 


ID 


Date 


Angular Separations (arcsec) 


Powerlaw Fitting for the western jet 


Eastern Jet 


Western jet 


Photon Index 


Flux (ergs cm 2 s - 1 ) 


1 


679 


2000 Jun 9 


21.5±0.5 








2 


1845 


2000 Aug 21 


22.8±0.5 








3 


1846 


2000 Sep 11 


23.4±0.5 








4 


3448 


2002 Mar 11 


28.6±0.5 


22.6±0.5 


1.75±0.11 


(1.9 ±0.4) x 10" 13 


5 


3672 


2002 Jun 19 




23.2±0.5 


1.71±0.15 


(1.6 ±0.3) x 10" 13 


G 


3807 


2002 Sep 24 


29.2±0.5 


23.4±0.5 


1.94±0.17 


(8.6 ± 1.5) x 10~ 14 


7 


4368 


2003 Jan 28 




23.7±0.5 


1.81±0.22 


(5.5 ± 1.0) x 10" 14 


8 


5190 


2003 Oct 23 




24.5±0.5 


1.97±0.20 


(3.1 ±0.6) x 10~ 14 
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4. JET MODEL 

4.1. Kinematic Model 

In the external shock model for afterglows of GRBs, the kinematic and radiation 
evolution could be understood as the interaction between the outburst ejecta 
and the surrounding ISM (Rees & Meszaros 1992). Microquasar jet systems 
are also expected to encounter such interactions. In this section, we describe 
our attempts after Wang et al. (2003) in constructing the kinetic and radiation 
model based on these models. 

We adopt the model of a collimated conical beam with a half opening angle 
9j expanding into the ambient medium with the number density n. The initial 
kinetic energy and Lorentz factor of the outflow material are Eq and Tq , re- 
spectively. Shocks should arise as the outflow moves on and heat the ISM, and 
its kinetic energy will turn into the internal energy of the medium gradually. 
Neglect the radiation loss, the energy conservation function writes (Huang, Dai, 
& Lu 1999): 

(r - l)M c 2 + a(T 2 sh - l)m sw c 2 = E , (1) 

where the first term on the left of the equation represents the kinematic energy 
of the ejecta, V is the Lorentz factor and M$ is the mass of the original ejecta. 
The second term represent the internal energy of the swept-up ISM, where T sh 
and 

^ SW 8-1*6 the corresponding Lorentz Factor and mass of the shocked ISM 
respectively, and m sw = (4/3)7ri? 3 m p n(0j/4). 

Coefficient a differs from 6/17 to 0.73 for ultra-relativistic and nonrelativis- 
tic jets (Blandford & McKee 1976). We adopt the approximation of a ~0.7 after 
Wang et al.(2003). Equation (1) and the relativistic kinematic equations 

dR P(T)c dR P(T)c 

{ dt )a 1 - f3(T)cos9 ,{ dt )r l + /3(r)cos# U 

can be solved and give the relation between the projected angular separation // 
and time t. In equations (2), the subscript 'a' and 'r' represent the approaching 
and receding jets in a pair of relativistic jets respectively. R is the distance 
between the jet and the source, which can be transformed into the proper motion 
separation by [i = R sin 9/5.3 kpc, and 9 is the jet inclination angle to the line 
of sight. We can get the [i — t curve numerically with the above equations. To 
be consistent with the work previously done to the eastern jet, we choose the 
same initial conditions that Tq = 3, Eq = 3.6 x 10 44 erg, and 9j = l.°5. Then 
the parameters needed to be fit are n and 9. 

In the case of the eastern jet, the number density of the ISM was assumed 
as a constant in the whole region outside the central source (Wang et al. 2003). 
This assumption does not work well in the case of its western counterpart. The 
western jet decelerated quite fast, requiring a local dense environment; however 
if the ISM is dense everywhere, the jet will be unable to travel that far from 
the central point source. As a result, we consider a model that the ISM density 
varies as the distance changes. For simplicity, we test the ideal case that the jet 
travelled first through a "cavity" with a constant velocity and then through a 
dense region where the jet was decelerated. A new parameter r, the outer radius 
of the cavity, is introduced and the ISM number density is set to be a constant 
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n outside this region and zero inside. The fittings improved a lot, but not well 
constrained because of the limited number of the data points. A combination of 
lightcurve fitting is required to further constrain the model parameters. 

4.2. Radiation Model 

In the standard GRB scenario, the afterglow emission is produced by the syn- 
chrotron radiation or inverse Compton emission of the accelerated electrons in 
the shock front of the jets (Wang et al. 2003 and references there). Wang et 
al.(2003) found that the reverse shock emission, originating from the electrons 
of the jet when a shock moves back through the ejecta, decay rather fast and 
describe the data of the eastern jet quite well. We thus take this model in our 
work as well. 

Assuming the distribution of the electrons obeys a power-law form, n^ c d^ c = 
Kj~ p d%, for 7 m < 7e < 7m , the volume emissivity at frequency v' in the 
comoving frame is given by (Rybicki & Lightman 1979) 

where F^v^v'J = tffrh F(x)x^/ 2 dx, with F(x) = x J +o ° K 5/3 (t) and 

K 5 /z(t) is the Bessel function. The physical quantities in these equations in- 
clude q and m e , the charge and mass of the electron, B±, the magnetic field 
strength perpendicular to the electron velocity, and v' m and is' M , the characteris- 
tic frequencies for electrons with 7 m and 7 M . 

Assuming the reverse shock heats the ejecta at time to at the radius i?o, the 
physical quantities in the adiabatically expanding ejecta with radius R will evolve 
as 7m = 7m(to)Ro/R,7m = 7m(to)Ro/R and K = K(t )(R/ R )-( 2 +p\ B ± = 
B±(to)(R/ Ro) 2 , where the initial values of these quantities are free parameters 
to be fitted in the calculation. With these assumptions, we can then calculate 
the predicted flux evolution of the jets. The comoving frequency v 1 relates to 
our observer frequency v by v = Du', where D is the Doppler factor and we 
have D a = 1/T(1 — (3cos9) and D T = 1/T(1 + j3cos9) for the approaching and 
receding jets respectively. Considering the geometry of the emission region, the 
observed X-ray flux in 0.3-8 keV band could be estimated by 

[^(-)ARD 3 j u/ ]d^ (4) 

where AR is the width of the shock region and is assumed to be AR = R/10 in 
the calculation. 

To reduce the number of free parameters, we set 7 m = 100 in our calculation 
because the results are quite insensitive to this value. We choose the time that 
the reverse shock takes place according to our kinematic model in section 3.1. 
Then we fit the data to find out the initial values of K and B± . 

Next step, we combine the kinematic and radiation fitting together. We 
know that the energy and the number density of the gas in the pre-shock 
and post-shock regions are connected by the jump conditions n' = C(T) n and 
e' = n(T)nm p c 2 , where C(T) and n(T) are coefficients related to the jet ve- 
locity. Therefore if we assume the shocked electrons and the magnetic field 
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Figure 2. Left: Proper motion fitting with asymmetric cavity. Parameters: 
9=68°, r c =14 arcsec, r w =18 arcsec, n c —0. 00675 cm~ 3 , n w =0.21 cm -3 . Right: 
Reverse shock emission fitting to the X-ray light curve of the jets. K c =0.09 
cm" 3 , # w =0.3 enr 3 , B c =0.5 mG, B w =0A mG. 



acquire constant fractions (e e and e B ) of the total shock energy, we have 7 m = 
e e {p - 2)m p (r - l)/(p - l)m e , K = (p — l)n' 7 r\ and B ± = VE^e 1 . 

If we further assume that the e c of the eastern and the western jets is the 
same, we may infer that K oc e 1 oc n for the two jets. As a result, we search for 
the combination of parameters that could satisfy the kinematic and radiation 
fitting, as well as the relationship K c /K Vi r^n c /n VJ . 

A set of parameters has finally been found (Please refer to the Left panel 
in Fig.2). The boundary of the cavity lies at r ~14" to the east and ~18" to 
the west of the central source. The corresponding number density of the ISM 
outside this boundary is ~0. 00675 cm -3 and ~0.21 cm~ 3 , respectively. Both 
values are lower than the canonical ISM value of ~1 cm' 3 , although the value 
in the western region is much higher than that in the eastern region. The electron 
energy fraction relationship is satisfied as K e /K w ^n c /n„ ~ 0.03. On the other 
hand, the other relation concerning the magnetic field strength could not be 
satisfied simultaneously by these parameters. Although the cavity radius and 
the number density are allowed to vary significantly, the best fitted magnetic 
field strength remains quite stable(~0.4-0.6 mG). One possible interpretation 
for this is that the equipartition parameter varies as the physical conditions of 
the jet varies; an alternative explanation may involve the in situ generation (or 
amplification) of the magnetic field. 



5. Conclusion and Discussions 

A GRB external shock model shows that a large scale cavity exists outside 
XTE J1550-564. This model has also been applied to another X-ray transient 
H 1743-322. Chandra X-ray and ATCA radio observations of H 1743-322 from 
2003 November to 2004 June revealed the presence of large-scale (~0.3 pc) jets 
with velocity v/c ~ 0.8 (Rupen et al. 2004; Corbel et al. 2005). Deceleration is 
also confirmed in this system. The external shock model describes the data of 
this source well. A cavity of size ~0.12 pc most likely exists, but the conclusion 
is not firm in this case. Even if there is no vacuum cavity, the ISM density is 
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found to be very low(~ 3 x 10~ 4 cm" 3 ), compared to the canonical Galactic 
value. 

These studies led us to the suggestion that in microquasars the interactions 
between the ejecta and the environmental gas play major roles in the jet evolu- 
tion and the low density of the environment is a necessary requirement for the 
jet to develop to a long distance. We find that microquasar jets can be classified 
into roughly three groups: small scale moving jets, large scale moving jets and 
large scale jet relics. For the first type, the "small jets" , only radio emissions are 
detected. The jets are always relatively close to the central source and dissipate 
very quickly, including GRS 1915+105 (Rodriguez & Mirabel 1999; Miller-Jones 
et al. 2007), GRO J1655-40 (Hjellming & Rupen 1995), and Cyg X-3 (Marti et 
al. 2001). The typical spatial scale is 0~0.05 pc and the time scale is several 
tenths of days. No obvious deceleration is observed before the jets become too 
faint. For the second type, the "large jets", both X-ray and radio detections are 
obtained, at a place far from the central source several years after the outburst. 
Examples are XTE J1550-564, H1743-322, and GX 339-4 (Gallo et al.2004). 
The typical jet travelling distance for this type is 0.2~0.5 pc from the central 
engine and deceleration is clearly observed. The last type, the "large relics", is 
a kind of diffuse structures observed in radio, optical and X-ray band, often ring 
or nebula shaped that are not moving at all. In this class, some well studied 
sources, Cygnus X-l (Gallo et al.2005), SS433 (Dubner el al.1998), Circinus X-l 
(Stewart et al. 1993), IE 1740.7-2942 (Mirabel et al. 1992) and GRS 1758-258 
(Rodriguez et al. 1992) are included. The typical scale for this kind is 1~30 pc, 
an order of magnitude larger than the second type. The estimated lifetime often 
exceeds one million years, indicating that they are related to previous outbursts. 

From these properties, it is reasonable to further suggest a consistent picture 
involving all the sources together. We make a conjecture that large scale cavities 
exist in all microquasar systems. The "small jets" observed right after the 
ejection are just travelling through these cavities. Since there are few or none 
interactions between the jets and the surrounding gas in this region, the jets 
travel without obvious deceleration. The emission mechanism is synchrotron 
radiation by particles accelerated in the initial outburst. The emissions of jets 
decay very quickly and are not detectable after several tenths of days. In some 
cases (e.g. XTE J1550-564), the cavity has a dense (compared to the cavity 
interior) boundary at some radius and the interactions between the jets and 
the boundary gas heat the particles again and thus make the jets detectable 
again. Those are the "large jets". The emission mechanism then is synchrotron 
radiation by the re-heated particles in the external shocks. Then, after these 
interactions, the jets lost most of their kinetic energy into the ISM gradually, 
causing the latter to expand to large scale structures, the "large relics", in a 
comparatively long time (several millions of years). 

The creation of the cavities is not clear at this stage. Possible mechanism 
may involve supernovae explosions, companion star winds or disk winds. Since 
some of the sources most likely have never had supernovae before and the winds 
from the companion stars are not strong enough, the accretion disk winds may 
be the most plausible possibility. However, these assumptions all require further 
observations to justify. 
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Microquasars are powerful probes of both the central engine and their sur- 
rounding environment. More studies of their jet behaviors may give us informa- 
tion on the ISM gas properties, as well as the ejecta components. It will provide 
insights of the jet formation process and offer another approach into black hole 
physics and accretion flow dynamics. 
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